Transformation optics with photonic band gap media 
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We introduce a class of optical media based on adiabatically modulated, dielectric-only, and 
potentially extremely low-loss, photonic crystals. The media we describe represent a generalization 
of the eikonal limit of transformation optics (TO). The foundation of the concept is the possibility 
to fit frequency isosurfaces in the /c-space of photonic crystals with elliptic surfaces, allowing them 
to mimic the dispersion relation of light in anisotropic effective media. Photonic crystal cloaks and 
other TO devices operating at visible wavelengths can be constructed from optically transparent 
substances like glasses, whose attenuation coefficient can be as small as 10 dB/km, suggesting the 
TO design methodology can be applied to the development of optical devices not limited by the 
losses inherent to metal-based, passive metamaterials. 
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WKB Approximation. 
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Transformation optics (TO) is a recently appreciated 
methodology for the design of novel electromagnetic de- 
vices pp. Devices created with the TO concept typically 
require exotic electromagnetic properties [2JG], which are 
unavailable in naturally occurring optical media such as 
solid-state crystals and amorphous substances (glasses). 
In fact, the required properties are so difficult to achieve, 
that TO-based designs would have been dismissed prior 
to the development of mesoscopic "metamaterial" fab- 
rication methods [2HI]. The desired material properties 
are usually available in metamaterials only at certain fre- 
quencies in the vicinity of a resonance band, which makes 
them inherently lossy. 

The lossy nature of metamaterials has impeded 
progress in certain areas of metamaterial and TO re- 
search, including sub-diffraction limited imaging (super- 
lensing) [5] and electromagnetic invisibility [TJ HJ EH • In 
the case of superlensing, which requires the exotic prop- 
erty of negative refractive index, an alternative, low- loss 
solution was proposed based on the existence of negative 
group velocity bands in photonic band gap media [7]- 
[9]. Photonic crystal-based imaging devices suffer from 
serious drawbacks, such as narrow operational band- 
width [10] and inefficient coupling between free space and 
the crystal [3J[9l[TT]. However, they offer the possibility 
of extremely low optical loss. So low the optical attenua- 
tion in a photonic band gap medium may be that it could 
span millions of wavelengths and yet be virtually trans- 
parent; the glass used in optical fibers [12] is a famous 
example of such a low- loss dielectric component. 

Transformation optics designs result in inherently com- 
plex media, typically requiring the exotic property of su- 
perluminal phase velocity (i.e., refractive index n < 1) 
combined with anisotropy (such as in beam shifters [TTJ 
[13]), spatial gradients in refractive index (as in the flat- 
tened Luneburg lens [14 ]), or both anisotropy and spa- 
tial gradients (as in the cloak of invisibility [US]). The 
"invisibility cloak" in particular has stimulated consider- 



able interest over the past several years, and serves as a 
natural test bed for various TO-inspired concepts and so- 
phisticated metamaterial designs. In this article, we sim- 
ilarly design a cloak to demonstrate the general ability 
to manipulate light propagation using photonic crystals, 
mimicking TO media. 

In its full generality, TO is able to map essentially any 
volume with an arbitrary bounding surface onto another 
arbitrarily shaped volume. For an invisibility cloak, one 
can restrict these transformations to radially symmetric 
maps. In two dimensions, a "perfect" cloak can be ob- 
tained from a coordinate transformation r' = q(r) that 
maps an annulus a < r < b onto the disk < r' < b [TJ [4] . 
To obtain a two-dimensional cloak, the constitutive pa- 
rameters in the physical space (with radial coordinate r) 
should be chosen as follows [4] [15] : 

e r = ^r = (q/r)/q, e e = fie = l/e r , ^z = u z = q(q/r\}) 

where q' = dq(r)/dr. The refractive indices correspond- 
ing to radial (n r ) and azimuthal (no) propagation in the 
perfect cloak (fij) are n r = q' and no = q/r, respectively, 
for both TE g] and TM [II [17] polarizations. 

In the geometrical optics limit one may abandon the 
notion of e and a in favor of the refractive index, since 
only the optical path length has relevance. The re- 
fractive index is defined as the wavenumber k corre- 
sponding to wave propagation in a certain direction, 
normalized to the vacuum wavenumber fco = cj/c, i.e. 
n(k) = |fc|/fco- Knowing the function n(k), or cj(fc)-the 
dispersion relation-is sufficient to describe the dynamics 
of light in transparent media [18]. Utilizing the disper- 
sion relation to describe wave propagation is more general 
than effective medium descriptions based on the tensors 
e, /i, since the properties of band gap optical media where 
propagating Bloch-Floquet modes exist, such as gratings 
and photonic crystals, can readily be analyzed [19] [20]. If 
the function n(k,r) changes slowly (adiabatically) with 



position 



r, i.e. 



Vf?n(fe, r) 



<C -r, light propagation can 




FIG. 1. (Color online) Dispersion relation for TM waves 
in a one- dimensional bi-layer photonic crystal with dielectric 
permittivities ei = 12, €2 = 1 and filling fraction t±/t = 
0.9. Horizontal axis is labeled by refractive index no = ke/ko 
corresponding to propagation in the periodic direction of the 
crystal, while the vertical axis is refractive index in the other 
direction (n r = k r /ko). 



be described in the WKB, or semiclassical, approxima- 
tion [H [21] . 

The adiabaticity of n(k, f) implies that it is continu- 
ous everywhere in space. Where a discontinuity exists, 
however, the function n(k,f) is not sufficient: it can- 
not describe transmission and reflection at the disconti- 
nuity, a limitation that applies to the sharp transitions 
between a photonic crystal and free space. Whereas an 
effective medium description would enable the reflectance 
and transmittance to be calculated from an effective wave 
impedance Z = y^e -1 , the wave impedance is generally 
an ill-defined concept for photonic crystals. The lack of a 
well-defined wave impedance in photonic crystals leads to 
the generation of multiple refracted and reflected beams 
(Bragg diffraction), and a dependence of reflectivity upon 
the microscopic details of the surface [5J [22]. Recent 
studies of negatively refracting photonic crystals with 
regards to their imaging properties show that improv- 
ing the coupling between free-space and crystal modes is 
possible [8] [9j [22] , although very difficult due to lack of a 
general theoretical model for the coupling. We note that 
the mode coupling problem is essentially unrelated with 
the dispersion engineering problem that we address here. 

Once the effective medium description is replaced by 
the uj(k,r) description, two features potentially useful 
for optical TO media can be noted. Firstly, materials 
with permittivity or permeability less than one, which 
are necessarily dispersive and lossy, are no longer needed. 
The dispersion relation of a photonic crystal near the T- 
point (k = 0) can be written as a Taylor expansion with 
even powers of k: 



uj 2 (k)=uj 2 +v 2 r k 2 r +v 2 k 2 + 0(k A ), 



(2) 



— ► 7,2 7,2 

or, equivalently, uj 2 {k)/c^ = ^ + -f + 0(k 4 ), where Cq 

2 _ 2 2 
is the speed of light in vacuum and v? r e = - ~" T -$p- . 



In the above we assumed that the propagation band 
has a positive group velocity, and thus, uj > cjp; sim- 



ilar formulas for a negative-refraction band can be ob- 
tained by changing the sign of v% G in the expressions 
above. The semi- axes of the elliptical contours of equal 
frequency in the /c-space are given by fcj? ax = n r uj/co and 
/-max _ neUj / C Q t i n the T-point (uj = cjp), no = n r = 0. 
Note that the effective indexes tlq^ are no longer related 
to any effective medium theory and, consequently, not 
restricted by Voigt-Reiss or Hashin-Shtrikman bounds 
arising from the classical homogenization theories [23] . 
Secondly, anisotropy of the refractive index is possible in 
dielectric-only crystals, and not only for TM waves, but 
also for TE waves. This contrasts with the behavior of 
TE waves in dielectric effective media, whose refractive 
index n = ^J~e~ z is always isotropic. 

While the approximately elliptic equal frequency con- 
tours ([2| are very common in the band structure of pho- 
tonic crystals in the vicinity of T-points (k = 0) [7] [TDJ 
[T9l [20] , some bands preserve their ellipticity throughout 
a relatively wide band. The existence of such an elliptic 
band in a one-dimensional photonic crystal formed from 
two alternating dielectric layers is illustrated by Fig. [I] 
The dispersion relation of the bi-layer crystal is known in 
analytical form [24]. Similar elliptic bands can be found 
in higher-dimensional photonic crystals, whose dispersion 
relation can be computed numerically. 

Theoretically, transformation optics is a powerful 
framework that enables unprecedented control over the 
propagation of electromagnetic waves. Anisotropy of 
the material properties implementing non-trivial coordi- 
nate transformations is a typical, and almost inevitable, 
property of the TO-based designs. A relatively simple 
TO-based device that uses anisotropic media is a beam 
shifter [13] . Its simplicity is due to the fact that a trans- 
verse shift of the beam can be obtained with a piecewise 
linear transformation, which results in piecewise constant 
material properties. In other words, the linear beam 
shifter can be implemented as a flat slab of a medium 
with uniform, although anisotropic, effective index. Pre- 
viously, reflectionless beam shifters based on embedded 
coordinate transformations have been demonstrated the- 
oretically in the effective medium regime [13] . 

A linear beam shifter with the shift parameter a = 
tan(</>) deflects the beam at an angle </>, resulting in 
a transverse shift Ay = ad, where d is the thickness 
of the shifter [13]. Below, we design and simulate a 
beam shifter operating in a higher Bloch-Floquet band 
of a photonic crystal. A bi-layer crystal illustrated 
and described in Fig. [I] is selected for this demonstra- 
tion. The operational frequency is chosen such that 
the two refractive indexes of the beam shifter medium, 

/ xl/2 

77,1,2 = f 1 + a 2 /2(l =p V(l + 4/a 2 ))j , match the two 
principal indexes of the photonic crystal; this leads to 
the choice k^t = 2.449, where t is the crystal periodic- 
ity. The resulting beam shifter has a shift parameter of 
a = 0.6187 and shift angle <\> — 31.7°. The principal val- 



ues of refractive index are rt\ =0.7374 and ri2 = 1.3561. 
The principal direction corresponding to smaller index 
n\ makes an angle 6 = 36° with the horizontal axis; the 
layers themselves make the same angle 9 with the vertical 
axis. 

The propagation of a Gaussian beam incident on a slab 
of homogeneous anisotropic medium with the refractive 
indexes described above is shown in Fig.[2Va). The beam 
traveling horizontally in free space propagates at an angle 
(j) = 31.7° inside the shifter. Note that the phase fronts 
remain vertical inside the shifter, indicating that the di- 
rections of the phase and group velocities differ. Fig.[2|b) 
shows the propagation of a beam with the same param- 
eters incident on a photonic crystal slab with the same 
effective indexes. The crystal slab is terminated with 
perfectly matched layers to reduce internal reflections. 
We observe that a free-space Gaussian mode excites two 
beams inside the crystal, indicating that its frequency lies 
within a higher propagation band. The upward-refracted 
beam follows the same path as the beam in Fig.[2Va). The 
downward-refracted beam is caused by Bragg diffraction; 
multiple solutions for the refracted wavenumbers can be 
understood as multiple intersections of a straight line 
k y = representing conservation of tangential wavenum- 
ber k y across the free space-crystal interface, with a peri- 
odically repeating pattern of isofrequency ellipses in the 
k x — k y space. To provide a clearer demonstration of the 
concept, we excite the crystal slab with a single mode 
of the crystal, modulated with the same Gaussian enve- 
lope, as shown in Fig. [2|c). The mode has no phase shift 
in the vertical direction, which corresponds to a horizon- 
tally incident beam. As expected, the second (diffracted) 
beam is eliminated while the first (refracted) beam is in- 
tact. Notably, the phase fronts of the Bloch modes in 
Fig. [2|c) remain vertical, in perfect agreement with ef- 
fective medium simulation in Fig. [2|a). 

To achieve a more complex TO design, such as a 
cloak, we can add an adiabatic spatial modulation to 
the one-dimensional crystal. The cylindrical cloak of in- 
visibility derived from equations (fil) requires no < n r . 
We have found numerically that in all elliptic bands of 
the bi-layer photonic crystal (except the lowest, acoustic 
branch), the refractive index for the propagation along 
the layers is always greater than in the other direction: 
n± < ni I ; this was found for both TM and TE polar- 
izations. Note that this property of higher-order bands 
is opposite to the behavior of the acoustic TM band, 
which has e±_ = (e -1 ) -1 < e\\ = (e) and therefore 
n±. — y^ > n \\ — v^L' as dictated by the Voigt-Reiss 
inequalities of homogenization theory [23 . This finding 
mandates that the crystal layers must be oriented ra- 
dially to achieve the required index distribution of the 
cloak, as shown in the inset to Fig. |3|a). 

The procedure for generating a radial transformation 
for eikonal cloaks has been introduced recently in 
Ref. [25 . In this case, it must be amended to account 




FIG. 2. (Color online) Parallel-slab beam shifter based on 
embedded coordinate transformations, (a) Gaussian beam in- 
cident from vacuum on a slab of anisotropic effective medium. 
Color shows out-of-plane magnetic field Re(H z ). (b) Gaus- 
sian beam incident from vacuum on a one-dimensional pho- 
tonic crystal. Color shows \H Z \. (c) Same as (b), except the 
Gaussian excitation is projected onto a single Bloch-Floquet 
mode of the photonic crystal. Color shows Re (H z ). 



for the additional geometric constraint. In contrast 
with effective medium-based cloaks, the unit cells 
in photonic crystal cloaks must change continuously 
between layers. In a rotationally-invariant design, to 
preserve the number of unit cells per circumference, 
unit cell size in the azimuthal direction, t, must change 
linearly with radius: t(r) — t^. To design a cloaking 
transformation, we assume that, besides t(r), there 
is one more parameter, p(r), that affects the band 
structure. Potential choices of p(r) in a bi-layer photonic 
crystal include (A) the dielectric contrast between the 
two components, ei/e2, (B) the filling fraction, ti/t, 
or (C) the scaling factor for both permittivities, a, 
such that ei,2(0 = e i,2(fr)c*(r), where a(b) = 1. To 
find the function p(r) that gives an eikonal cloak, one 
may start from the differential equation of Ref. [25] , 
n r (p,t)-n 9 (p,t) i wnere n r,e are known functions of 



dr 



p and t calculated from the dispersion relation. As- 
suming a linear function for t(r), the unknown function 
p(r) is found by integrating a first-order ODE ^ = 

/ n r (p,t(r))-ng(p,t(r)) dn e (p,t(r)) dt\ ( dn e (p,t(r)) \ ~ 1 ^ 

y r dt dr J \ dp J 

F(p,r). The choice (C) allows an analytical simplifica- 
tion of the dispersion engineering problem. When both 
the dielectric contrast and the filling ratio are constant, 
the band structure is a function of a single parameter, 
(r) = t(r)a x / 2 (r), and the transformation can be found 



as a simple quadrature: ff = J nr{ ^ neW ff- The 
resulting cloaking transformation and permittivity dis- 
tribution for photonic crystals with radius-independent 
dielectric contrast 61/62 = 12 and several choices of 
the filling ratio are plotted in Fig. |3ja). In the re- 
mainder of this article, we will use the filling fraction 
t\jt = 0.9, which leads to an eikonal cloak with aspect 
ratio a/b = 0.2631. Fig. [3lb) shows the distribution of 
refractive indexes in such a cloak. 
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FIG. 3. (Color online) (a) Scaling of the dielectric constants 
with radius, for the photonic crystal cloaks with ei/e2 = 12 
(TM polarization). The inset shows the geometry of the cloak. 
(b) Refractive index distribution in the photonic crystal cloak 
with dielectric contrast ei/e2 = 12 and filling fraction t±/t = 
0.9. 




FIG. 4. (Color online) (a) Cloaking of a normally incident 
Gaussian beam modulated with the desired mode profile of 
the photonic crystal, (b) Same as (a), but with the angle of 
incidence 0i nc = 15°. 



The operation of the photonic crystal cloak is demon- 
strated using electromagnetic simulations performed 
with COMSOL, a Finite Element Method solver. The 
geometry of the cloak consists of two concentric circles of 
radii b = 0.4 (m) and a = 0.10524 (m). The inner disk 
represents the cloaked object, on which perfect electric 
conductor boundary conditions are imposed. The space 
between the circles is filled with a quasi-periodic photonic 
crystal with N — 512 radially converging layers. In order 
to suppress additional beams, a single mode of the crystal 
is injected in the simulation in the same fashion as was 
done in the beam shifter model in Fig. [2|c). As seen in 
Fig. Qa), a Gaussian-like beam is created in the crystal. 
The beam begins to diverge as it travels towards the cen- 
ter; then it bends around the cloaked object and re-forms 
as another beam on the opposite side of the cloaking de- 
vice. Excitation of the same mode with a non-zero phase 
shift per unit cell results in a Gaussian beam traveling at 
an angle, as shown in Fig. Kb). 

In conclusion, we have demonstrated theoretically and 
in high-resolution electromagnetic simulations that the 
concepts of transformation optics can be implemented 
with Bloch-Floquet waves in adiabatically modulated 
photonic crystals. This includes a demonstration of a 
dielectric-only optical cloak, whose practical implemen- 
tation can be composed of materials that are macroscop- 
ically transparent in the entire visible spectrum. This 



paradigm creates hope for a new class of optical devices 
capable of being scaled to arbitrarily large dimensions. 
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